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Abstract 

Using a computer model for the crystal extraction, we inves- 
tigate the bent-crystal parameters optimal for the extraction ex- 
periment. The optimal crystal curvature is found to be near 1 
GeV/cm (for pv/R), i.e. a factor of 2-3 higher than for the 
crystal application in beam lines. An influence of the accelerator 
optics on extraction is discussed. A possibility of using the high-Z 
crystals for extraction is considered. The simulations for the on- 
going experiments at the CERN-SPS and the Fermilab Tevatron, 
and for the proposed extraction at LHC, are presented. 
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1 Optimization of single-pass channeling 



Optimization of channeling in a single pass is well understood]]]]. Two effects 
cause the particle dechanneling in a bent crystal: the bending dechanneling 
due to centrifugal force, and the scattering. For a given bending angle Ob, 
a larger crystal length L means smaller curvature 1/R=Qb/L and hence 
smaller bending dechanneling. However, the L increase means larger loss 
for scattering (characterized by the 1/e-dechanneling length Ljj). There is 
an optimal L, when losses from bending and scattering are roughly equal. 
In a harmonic approximation for the interplanar potential, the bent-crystal 
transmission can be estimated as 

We used the notation p=R c /R (R c is the critical radius). As is the number 
of particles trapped in a straight crystal. The factor (1 — R c /R) 2 describes 
the As reduction due to bending dechanneling. The exponent describes a 
reduction due to scattering; Oe>=Ld/R c is a constant defined by the crystal 
properties. In a bent crystal, Lp is reduced by (1 — R c /R) 2 , just the same 
factor as for bending dechanneling. See review [I]] for further details and 
references. The dislocation dechanneling in crystals is reviewed in Ref. ||. 

2 Optimization of multi-pass channeling 

2.1 Qualitative discussion 

If a crystal works in the extraction mode at a circular accelerator, the cir- 
culating particles may pass through it many times. The simplest estimate 
of the total (with multipasses) efficiency of the particle extraction F may be 
derived using the following two simplifications: 

(a) the probability A\~ of the particle capture into the channeling mode 
in the fe-th passage through the crystal is the same as that probability in 
the first passage, Ak = A\ = A (actually A& decreases with k as the beam 
divergence grows with the scattering in crystal); 

(b) the probability qk of the particle loss, either in the crystal or at the 
vacuum chamber wall, with the fc-th passage through the crystal is the same 
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as that probability in the first passage, qk = qi = q (actually qk may increase 
with k if the particles reach the vacuum chamber wall with scattering). 

The circulating particles are gradually removed from the ring with these 
two processes - crystal channeling and nuclear interaction, with respective 
probabilities A and q per a single passage of crystal. The number of the 
particles remaining in the ring after k passages through the crystal is Nk = 
N (l— A) k (l— q) k if the initial number was N Q . The particle number extracted 
by the (k + l)-th passage equals F k N Q = AN {1 - A) k {l - q) k , and total 
efficiency amounts to 

F = Y F k = : : exp(-L/L D ) (2) 

^ k A- (1 -q) + q K J 

The exponent takes into account dechanneling of the trapped particles. Since 
both assumptions (a,b) tend to overestimate F, Eq. (0) is the upper limit for 
the real efficiency. 

Eq. (0) well agrees, e.g., with the computer simulation[Q of proton ex- 
traction from a 20-TeV SuperCollider, Fig. 1. The analytical curve (|j) was 
predicted in Ref. || for 90-/zrad bending; the points are from the more-recent 
simulation for 100-/irad bending^. 

To understand Eq. (H) better, it is useful to consider several limiting cases. 

1. If the probability q of the unchanneled particle interaction with nuclei 
over the crystal length is very high, q ~1, the secondary passes are 
not possible. Eq. (|J) is reduced then to the efficiency of a single pass: 
F ~ Aexp(-L/L D ). 

2. Typically, q ~ L/L^ —0.1 over a single pass in crystal, where 

is the nuclear interaction length. If the probability of capture into 
channeling mode is very small, A <C q, Eq. (|2]) is reduced to F ~ 
(l/g)Aexp(— L/Ld). This is simply the single-pass efficiency times the 
mean number of passes (1/q) before a nuclear interaction. 

3. The opposite case q A would mean that the circulating particles are 
removed from the ring mostly with the crystal channeling rather than 
with nuclear interactions. The efficiency is then close to 1, with Eq.(Q) 
reduced to F~ (1 — q/A) exp(—L/L D ). 

A shorter crystal disturbs the particles less. The amplitude of betatron 
oscillation of a scattered proton grows with the scattering angle squared 
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Figure 1: Extraction efficiency for a 20-TeV collider as a function of the 
Si(110) crystal length. The curve (•) is analytical prediction Eq. (2). The 
other data are from computer simulation [4]: the first-pass ((g)) and the overall 
(o) efficiency 

9 2 '- x max=(xo + P 2 ^) 1 ^ 2 , P i s the accelerator function, xq the crystal edge 
position. When x max exceeds the aperture, the proton is lost. Keeping 9 2 s 
below some critical value, one may reduce L to get a respective increase in 
the allowed number of passes ~1/L before a loss. Furthermore, the scattering 
over L is reduced, so the scattered protons have smaller divergence in further 
passes. 

On the other hand, a short crystal need to be more strongly bent to keep 
@£j. With pv/R greater than ~1 GeV/cm in silicon, the crystal acceptance 
rapidly decreases. This cuts the extraction efficiency at too-small L. 

2.2 Optimization of crystal parameters 

Detailed computer study of extraction should take into account the accel- 
erator optics (betatron motion, aperture) and the processes of channeling, 
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scattering, and nuclear interaction in crystal. These processes depend essen- 
tially on L. 
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Figure 2: The SPS extraction efficiency simulated as a function of the 
Si(110) crystal length. For perfect surface (o) and t—1 /im (•). The dots 
(<8>) are for the U-shaped design and t=20 /im. Also shown is the measured 
range of efficiencies, 10-15 % for the 4-cm U-shaped crystal. 



Figs. 2-4 show the extraction efficiency simulated as a function of L, for 
the experiments at 120-GeV SPS |, |, 0, §, 900-GeV Tevatron g, [Eg, and 7- 
TeV Large Hadron Collider (see refs. in ||). The accelerator setting and the 
computer model have been earlier described in Refs. [11], [12], |13| respectively. 



In Figs. 2-3 two cases were studied: an ideal crystal, and the crystal with 
imperfect surface. In the latter case an inefficient layer of some thickness 
t appears near the crystal surface, defined by the surface irregularities. In 
Fig. 4 we consider imperfect surface only. Also shown in the SPS and LHC 
cases is the simulation for a realistic model of the "U-shaped crystal" || [7|, § 
which included the end straight parts (5 mm for L>5 cm and 1-4 mm for 
L<4 cm), and scattering in "the legs of U" and in the bending device. The 
incident beam was nearly parallel at the SPS and LHC, but divergent (11.5 
/xrad rms) at the Tevatron. The efficiency measurements have been reported 
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so far from the SPS for the crystals of the 3-4 cm length only. The simulation 
for the U-shaped crystal 4-cm long in realistic model has given the efficiency 
in good agreement with the experimental values of 10-15 %, Fig. 2. Notice 
that an agreement was found only in the assumption of imperfect edge of 
crystal. This edge imperfection (the first-pass suppression) predicted |11 
for the SPS extraction was confirmed by the recent experiment ||. Other 
comparisons of the expected and measured parameters in crystal extraction 
see in Refs. |, g, 0. 




Figure 3: The Tevatron extraction efficiency for the ideal (o) and imperfect 
(•), t=l fim, crystals as a function of the Si(110) crystal length. 



In all the Figs. 1-4, the optimum is found near pv/R ~1 GeV/cm. The 
short-crystal optimum is most spectacular at lower energies, ~120 GeV, 
where scattering is appreciable. 

All the circulating protons with sufficient amplitudes are intercepted by 
the crystal sooner or later, irrespective of its transverse size. Therefore this 
parameter is not so important for extraction. 
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Figure 4: Extraction efficiency for 7-TeV LHC with Si(llO) (o) and Ge(llO) 
(•) "U-shaped" crystals (imperfect surface). 

2.3 Extraction with high— Z crystals 

Crystals with higher Z than in Si (Ge and possibly W) have evident benefits, 
such as higher critical angle O^Z 1 ^, and smaller R c ~l/Z; <d D ^Z 2 ^. There- 
fore the application of high- Z perfect crystals for a "single-pass" channeling is 
advantageous. Interestingly, in crystal extraction (which is multi-pass chan- 
neling essentially) this may not be true. The primary incident particles can 
be practically parallel due to very small impact parameters b (in the range 
of 1 ^m) at crystal; then 9 C doesn't matter. In multipass mode, one has 
to trap the protons scattered in the crystal. Then the scattering angle and 
interaction length do matter as well as the channeling properties. 

In general, one could consider a constant crystal lattice and vary Z of the 
atoms (we keep Z/A^const below). Then 9 C is scaled as Z 1 ' 3 . The scattering 
angle per unit length is scaled as Z. Increasing Z, one scales the trapping 
efficiency as Z~ 2 I Z (i.e., the ratio of 9 C to particle divergence). Since R c ~l/Z, 
one can shorten L~l/Z. The scattering angle per a single passage of crystal 
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is ~Z 1 / 2 , and the efficiency is Z^ 1 ^ (decreasing with Z). 

The final answer is more complicated, as one takes into account the 
dechanneling loss (reduced in shorter crystal: L/L^ ~ Z~ 2 / 3 ), and absorp- 
tion length L N r^Z~ 2 ^ 3 (absorption loss is also reduced: L/L N ~ Z -1 / 3 ). 
Fig. 4 shows the simulated efficiency of proton extraction from LHC by 
Ge(llO) crystal in the same conditions as for Si(llO). As expected, the 
efficiency with Ge crystal has slightly decreased compared to Si. Therefore, 
high-Z crystals are hardly advantageous in extraction. 



2.4 Accelerator parameters 

The channeling efficiency is defined by the divergence of incident particles 
w.r.t the crystal planes. In the "diffusion mode" of extraction, protons are 
intercepted by crystal with very small impact parameters b in the first pass. 
Hand-in-hand, the primary divergence in the "diffusion plane" is also small: 
b' <C 9 C . The crystal may deflect protons either in this plane (Protvino- 
CERN case [14, ||) or normally (Fermilab[[|). In the latter case, primary 



divergence is scaled as (3~ 1 / 2 . Interestingly, in the Protvino-CERN case the 
efficiency hardly depends on j3. The starting divergence of protons b' <C 6 C ] 
the secondary divergence is defined by the scattering in crystal, which is 
independent of the machine optics. 

Increasing (3 one could increase primary b to avoid a "septum width" t 
due to crystal surface irregularities (and get an efficient first-pass channeling). 
However, recent measurement s fl"5|| indicate t~50 /zm, while b is well below 1 
/im. In that case, the extraction becomes insensitive to (3. 

The accelerator a-function sets a correlation of the mean incident angle 
x' of particles with the transverse coordinate: x' mean =— ax/ (3. The angle 
of incidence varies by Ax'=at cr / (3 over the crystal thickness t cr . E.g., with 
t cr =3.5 mm in the SPS case, the x' variation over the crystal face becomes 
2.07x3.5mm /90nr»80 /irad > 6> c =14 /zrad! A thick crystal (t cr > (39 c /a) 
should loose in efficiency. 

We leave aside the discussion of the betatron tune effects in crystal extrac- 
tion. Two points can however be mentioned. At the rational values of tune 
the trajectory of a circulating particle becomes periodical; if the first pass in 
crystal was through an inefficient layer, some later passes suffer the same way 



resulting in a dip of efficiency at the rational tunes [11|. Note however that 



in a real machine the other strong effects are essential at the rational tunes. 
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Another relevant point is that crystal theorists predict [IB, [T3| that in an 
aligned bent crystal the unchanneled particles are deflected in average at the 
angles of order —9 C , opposite to the crystal bending. In terms of the betatron 
motion, a deflection at 8 C is a betatron phase advance of AQ fh (38 c /xq per 
each traverse of crystal. Typical AQ ~0. 1-0.2 is big and should be taken 
into account in a study of the tune effects in crystal extraction. 



2.5 Pre-scattering 



To raise primary b, some ideas were proposed (T7[ and employed ]Tjfl to pre- 



scatter protons in a thin target before they hit a crystal. In a target of length 
s, the scattering angle is 9 S =(E S / E)(s/ Lr) 1 / 2 ; E s =14 MeV, Lr radiation 
length. Then, the impact parameters in the first pass are h ~ /3 2 8 2 /2xq |L^] , 
which gives a relation for s: 

2bx L R E 2 

s - wm {) 

Notice, s ~ E 2 . E.g., requiring b >50 jum at crystal, one needs a Si target as 
long as s>l mm at 120-GeV SPS, but s>40 mm at 900-GeV Tevatron. For 
comparison, an effective length of the first pass at SPS is ~5 mm (reduced 
near the bent face); in the 900-GeV case, the 40 mm is as long as the crystal. 
Therefore, an amorphous pre-scatter might work at ~100 GeV and below. 
At >1 TeV this idea can no longer be applied. Notice that a pre-scatterer 
must be positioned (closer to the beam than crystal) with accuracy ~6 (~50 
/im)! Otherwise, protons traverse it many times, and the benefits are lost. 
The pre-scatter becomes more interesting with a crystalline scatterer 



producing angular kicks ~# c . The minimal length required is A/4, with A 
being the channeled-particle wavelength (A~0.1 mm at 1 TeV in Si). The 
drawback is the need to tune two extra parameters: position (with accuracy 
~10 yum) and angle (with accuracy <9 C ) of the pre-scatterer. 

The alternative to the above ideas is the intrinsic scattering in the edge 
of bent crystal. It was shown with simulations for 120-GeV, 900-GeV and 
7- TeV extraction experiments that a difference in efficiency between an ideal 
crystal and a crystal with imperfect surface vanishes with increasing energy 
i?jn], |T2], Therefore, the problem may disappear in TeV range. 
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3 Conclusion 



The computer study of the crystal-extraction at accelerators indicates an 
interesting physics of "multipass channeling". The optimal extraction is 
achieved with a short crystal, and is due to the multipass mode when the re- 
duction in a single-pass transmission is overcompensated with the increased 
mean number of passes. The multipass mode may well be the only feasible 
one at high energy accelerators, due to infinitesimal b. This mode must be 
studied (at the optimum suggested) as it provides a deeper insight into the 
crystal-extraction physics. Different ideas for studying the multipass chan- 
neling at circular accelerators have been reviewed in Refs. ]TT] , |T2"| . 
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